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Abstract. We present the first results of the treatment of grain growth in our 3D, two-fluid (gas+dust) 
SPH code describing protoplanetary disks. We implement a scheme able to reproduce the variation of grain 
sizes caused by a variety of physical processes and test it with the analytical expression of grain growth 
given by Stepinski & Valageas (1997) in simulations of a typical T Tauri disk around a one solar mass star. 
The results are in agreement with a turbulent growing process and validate the method. We are now able 
to simulate the grain growth process in a protoplanetary disk given by a more realistic physical description, 
currently under development. We discuss the implications of the combined effect of grain growth and dust 
vertical settling and radial migration on subsequent planetesimal formation. 

1 Introduction 

Planets are thought to form in disks around young stars through the growth of dust grains (Dominik et al. 2007) . 
Collisions and aggregation govern the first steps from micron-sized particles to decimetric pre-planetesimals, 
which then form kilometric planetesimals via still-discussed mechanisms, and ultimately planets. Observations 
of protoplanetary disks support this mechanism by showing evidence of dust grain growth (e.g. Apai et al. 
2004). The understanding of this physical process is necessary to describe the first steps of planet formation. 
We describe here the implementation of grain growth in our simulations of protoplanetary disks. 

2 Physics of protoplanetary disks 

Protoplanetary disks are mostly composed of gas and typically contain about 1% of their mass in dust. The gas 
evolves in a pseudo keplerian rotation due to its radial pressure gradient, which decreases the effect of the central 
star's gravity. The energy dissipation due to the gas viscosity induces an inwards mass flow and an outwards 
angular momentum flow. Dust grains only feel the gravity field of the central star, but are not affected by either 
pressure gradient or viscosity and, if alone, would therefore be in keplerian rotation. 

However, gas and dust interact via a drag force which depends on their relative velocity and slows down 
the dust, causing it to migrate radially inwards and to settle vertically towards the midplanc. Those processes 
are strongly dependent on the size s of the grains: the largest particles are totally decoupled from the gas, the 
smallest ones are so strongly coupled that they follow the gas motion, and in the intermediate regime, particles 
undergo settling and migration with a varying efficiency (Barriere- Fouchet et al. 2005). 

Finally, the size of dust grains can vary: depending on their relative velocities and material properties 
(Dominik & Tielens 1997), collisions between solid particles can make them stick and grow, or conversely break 
them into smaller pieces. 
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3 Grain growth in protoplanetary disks 

A first description of the grain growth process is given by Stepinski & Valageas (1997). They model a turbu- 
lent, vertically isothermal protoplanetary disk, in which gas and dust are represented by two separate phases 
interacting via aerodynamic drag in the Epstein regime. Their solid particles are supposed to stick perfectly 
during collisions, and can therefore only grow. The variation of their size s is given by the following analytical 
expression: 

^ = v/2^o^ a^EI, (3.i) 

dr pd be 

where Ro is the Rossby number for turbulent motions, a the Shakura & Sunyaev (1973) viscosity parameter, 
p s the density of matter concentrated into solid particles, p,i the intrinsic density of the grains, C s the local gas 
sound speed, and Sc the Schmidt number of the flow which estimates the effect of gas turbulence on the grains. 
Sc is defined by 

Sc= (l + r» k < s )^l + ^, (3.2) 

where fik is the local keplerian velocity, t s the dust stopping time, v the mean relative velocity between gas and 

ds 

dust, and V t a turbulent velocity. The growth rate — depends on s through the stopping time 

at 

t s = ^r, (3.3) 

Pg O s 

where p g is the gas density. 

The motion of a dust grain, as mentioned in Sect. [21 and the grain growth process itself presented here, 
which both depend on the grain size s, are coupled phenomena. The understanding of the global evolution of 
dust in disks therefore requires a numerical treatment. 



4 Grain growth with an SPH code 

Our 3D, bi-fluid, Smooth Particle Hydrodynamics (SPH) code has been developed over recent years to model 
vertically isothermal, non self-gravitating protoplanetary disks. The effects of radiative transfer, magnetism, 
self-gravity and turbulence are neglected. Gas and dust are treated as two separate phases and are coupled by 
aerodynamic drag. The code and its results on dust migration and settling are presented in Barriere-Fouchet 
et al. (2005). It has also been applied to grain stratification in GG Tau's circumbinary ring (Pinte et al. 2007) 
and to gaps opened by planets in the dust phase of protoplanetary disks (Maddison et al. 2007; Fouchet et al. 
2007). 

The assumptions we have made are very similar to those of Stepinski & Valageas (1997) mentioned in Sect.[3j 
their prescription for grain growth is therefore easy to implement in our code. In contrast to our previous work, 
we now allow the grain size s assigned to each SPH particle, assuming it to represent the typical size of dust 
grains at its position, to vary with time following Eq. 13.11 We take the initial grain size distribution to be 
uniform. 



5 Results 

We model grain growth in a typical T Tauri disk of mass Mdisk = 0.01 M®, with a total dust mass Mdust — 
0.01 Mdisk, extending between 3 and 300 AU from the central star of mass M* = 1 Mq. We ran simulations with 
100,000 SPH particles, evolved over 40,000 years, for a series of initial dust grain sizes ranging from sq = 1 pm 
to 1 mm. Figures [T] and [2] show the size distribution in the disk as a function of the radial distance to the central 
star at the end of the simulation with s = 1 A* 111 - 

We find that grain growth occurs very quickly, especially in the inner disk where the density is the highest 
(on the order of 10 -10 kg.m -3 ), as one would expect from Eq. (|3.ip . In this region, the dust grains reach 
centimetric size (see Fig. [IJ in only a few timesteps. In the outer parts of the disk, where the density is far 
lower (as low as 10 -15 kg.m -3 ), the grains grow much more slowly and their size stays below the millimeter 
(see Fig. 
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Fig. 1. Radial dust grain size distribution at the end of the SPH simulation with so = 10 /jm, in a linear scale. 



At any place in the disk, the growth timescale is far smaller than the migration timescale. Hence, after 
a transitory regime where grains grow from their uniform initial size, a quasi-stationary regime quickly takes 
place. Indeed, an inwards migrating grain coming into a denser region grows almost instantaneously to a size 
characteristic of the grains at its new position, which in turn evolves slowly over time. The global shape of the 
plots shown in Figs. Q] and stays the same as time goes on, but the size distribution slowly progresses to larger 
sizes. 

The initial size so has very little influence on the final size distribution, and in particular on the maximum 
size reached. Indeed, for a given dust disk mass, the smaller the dust grains, the more numerous they are, and 
the more frequent their collisions are. Consequently, even if the initial grains are very small, they will very 
quicky reach the asymptotic regime and lead to centrimetric grains in the central parts. 

Our results are consistent with those obtained by Dullemond & Dominik (2005). With a model solving the 
coagulation equation in presence of turbulence, they found a very fast grain growth in T Tauri protoplanetary 
disks, depleting very small sizes and producing centimetric grains. 



6 Tests with a semi-analytical model 



In the Stcpinski & Valageas (1997) prescription, presented in Eqs. (|3.1[) and (|3.2j) . the growth rate depends on 

the local differential velocity between gas and dust. This quantity is determined by the SPH code. In order 

to test the validity of our implementation, we develop a semi-analytical model of viscous disks, to which we 

compare our SPH simulations. In this model, as in our code, the disk is assumed to be vertically isothermal, non 

turbulent, non sclf-graviting and all magnetic effects arc neglected. To describe the gas migration, the viscous 

character of the gas has to be considered. For simplicity, the gas flow is incompressible and in a stationary 

regime, and the migration timescale is assumed to be large compared to all other timescales. Under these 

v 

hypotheses, the gas radial velocity is proportional to , where v is the gas kinematic viscosity. 
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Fig. 2. Same as Fig[T]m a logarithmic scale. 



The equation of dynamics for the grains is: 

dvd 
dt 



Vd 



(6.1) 



where the d and g subscripts refer to the dust grains and the gas, respectively. For small grains, we derive an 
expression for the components of the differential velocity (Laibe & Gonzalez 2008). Taking into account the 
initial conditions for the dust, which is injected with the same local velocity than the gas, i.e. Vd(t = 0) = 
v g (t = 0), we obtain 





V g,r 




Vd.e 




= 


Vd,z 


- %,* 


= 



r 




,-t/t. 



(6.2) 



v 

Using the — dependence of the gas radial velocity, we have rewritten the viscosity term so that it appears only 
r 

v 2 r 

implicitly in the term. Therefore, both the viscosity and the radial pressure gradient contribute to the 
r 

radial differential velocity. 

We ran a new simulation with the same settings as in Sect. [5l but this time v in Eq. f|3. 2|) is computed 
via Eq. (|6.2p instead of using the individual gas and dust velocities. Figure [3] displays the resulting grain size 
distribution, it resembles the one shown in Fig.[Tl in particular the maximum grain size reached at the inner disk 
edge is of the same order of magnitude, about 6 cm. This provides a consistency check of our implementation 
in the SPH code. The residual difference between the two simulations is due to the restricting hypotheses made 
in the analytical model. 

In order to investigate the relative importance of the two terms contributing to the radial differential velocity, 
we ran a simulation omitting the viscosity term. As illustrated in Fig. |4l the maximum grain size in the inner 
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Fig. 3. Radial dust grain size distribution obtained by the analytical model. 



disk is about twice as large as the one obtained in both previous simulations. This shows the importance of the 
viscosity term in Eq. (|6.2p . and that it should not be neglected. 

7 Conclusion 

The physics of grain growth in protoplanetary disks is complex and involves coupled processes, it therefore 
requires the use of numerical simulations. We have implemented in our 3D SPH code a mechanism able to 
treat this grain growth and validated it through the use of the simple model of Stepinski & Valageas (1997) and 
comparisons with a semi-analytical model we developed. 

In accordance with physical intuition, dust grain grow much more quickly in the denser, central regions of 
the disk, where ccntrimctric sized are reached. The growing time is smaller than the migration time and a 
quasi-stationnary distribution of grain size appears in the disk. The small grains are depleted too rapidly to 
be consistent with observations of protoplanetary disks, showing the need to take into account other processes 
such as shattering of large grains in high-velocity collisions. 

In order to treat a more realistic grain growth process, we are working on a more detailed model taking into 
account microscopic interactions between the grains, kinetic energy dissipation, porosity, and re-fragmentation, 
from which we will derive a prescription to include in our code in the same way as the one we have presented 
here. We also plan to add a simple, but realistic, treatment of turbulence. 
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Fig. 4. Radial dust grain size distribution obtained by the analytical model when the viscosity term is neglected. 
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